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Positive definiteness of the quadratic part of the action of the Hawking- Turok instanton is 
investigated. The Euclidean quadratic action for scalar perturbations is expressed in terms of a 
single gauge invariant quantity q. The mode functions satisfy a Schrodinger type equation with a 
potential U. It is shown that the potential U tends to a positive constant at the regular end of 
the instanton. The detailed shape of U depends on the initial data of the instanton, on parameters 
of the background scalar field potential V and on a positive integer, p, labeling different spherical 
harmonics. For certain well behaved scalar field potentials it is proven analytically that for p > 1 
quadratic action is non-negative. For the lowest p = 1 (homogeneous) harmonic numerical solution of 
the Schrodinger equation for different scalar field potentials V and different initial data show that in 
some cases the potential U is negative in the intermediate region. We investigated the monotonously 
growing potentials and a potential with a false vacuum. For the monotonous potentials no negative 
modes are found about the Hawking- Turok instanton. For a potential with the false vacuum the 
HT instanton is shown to have a negative mode for certain initial data. 



PACS: 98.80.Hw, 98.80.Cq ;r-qc/9804056 BUTP-98/12 



I. INTRODUCTION 

A recent proposal (!]] by Hawking and Turok (HT) to use singular instantons to describe creation of an open 
inflationary universe caused a cascade of works (see e.g. Originally the HT instanton was discussed in a model 

with a scalar field coupled to gravity. Later it was argued that the inclusion of the four form field could lead to a 
solution of the cosmological constant problem || . 

It is assumed that the quadratic action of the HT instanton is positive definite On the other hand it is 
known that some Euclidean solutions such as the bounce in a flat space-time, the Hawking-Moss instanton, and the 
Giddings-Strominger wormhole have the negative modes Jl0|-p2[ . 

The aim of the present note is to investigate positive definiteness of the quadratic action for the HT instanton. We 
shall analyze the Euclidean quadratic action which is obtained by the analytic continuation of the Lorentzian quadratic 
action for a scalar field in a closed universe recently derived in [fl3|| . After the Hamiltonian reduction, the quadratic 
action can be expressed in terms of a single gauge invariant quantity q. Expanding the perturbations by harmonics on 
the unit 3— sphere the quadratic action can be represented as a sum of quadratic actions for different harmonics. The 
mode functions satisfy a Schrodinger type equation, with the potential U depending on the background (instanton 
solution) and harmonics number p. We shall demonstrate that the potential tends to a positive constant at one 
(regular) end of the instanton. The detailed shape of U depends on the instanton initial data, on parameters of the 
background scalar field potential, and on the harmonic number p. For p > 1 we prove analytically that quadratic 
action is non-negative. General proof is not working for the lowest p = 1 (homogeneous) harmonic. A numerical 
solution of the Schrodinger equation shows that in some cases potential U is negative in the intermediate region. For 
the monotonous (quadratic and quartic) scalar field potentials we found that the zero energy wave function has no 
nodes. We investigated a case of a potential (with a false vacuum) in which the Coleman-De Luccia (CD) solution 
ful exist. We found that the HT instanton for a certain initial data (which correspond to "overshooting" the CD 
solution) has a negative mode. 

The rest of the note is organized as follows. In the next section we discuss main properties of the HT instanton. 
In section 3 we consider the Euclidean quadratic action. In section 4 we present our numerical results. Section 5 
contains the concluding remarks. 
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II. THE HAWKING-TUROK SINGULAR INSTANTON 



The Euclidean action for a scalar field coupled to gravity has the form 

S E = J(-^R + l(d^) 2 + V( v ))^d 4 x, (1) 



where k — 8irG is the reduced Newton's constant and V(tp) is a scalar field potential. 
We shall be interested in 0(4) invariant solutions and metric is parameterized as 



ds 2 = da 2 + b 2 (a)dnl , (2) 



where dft 2 = dip 2 + sin 2 (%li)dVt\ is the line element of the unit 3-sphere. 
The field equations are 

h 1 r)V k 

^" + V' = ^' i%' 2 + ^)), (3) 

where prime denotes d/da. 

The HT instanton has the following initial conditions at a — ► |jj 

<p(0)=<f>o, ¥>'(0) = 0, &'(0) = 1, (4) 

which implies 

1 f)V K 
<p{cx)= Vo + ^+0{a% b(*)=a--V& Q )<J 3 + 0(a 5 ), (5) 

where ipo is some constant. 

If one starts at a = with the initial conditions Eq. (^), the scale factor b grows, reaches maximal value and 
decreases again developing the second zero at some cr = er/. Next HT jlj argue that as we approach <jf the potential 
terms become irrelevant in the field equations (JsJ) : the first equation then implies that (p 1 oc b~ 3 and the second yields 
b OC ((7/ - cr) 1/3 or 

V = Vf-\jl^H°f-°), b = C(a f -a) 1 / 3 , (6) 
V o/s^ 

where ipf and C are some constants. This assumption is certainly true for potentials, which grow sufficiently slowly 
as tp — > oo, namely 

V( V (<7/ - a)) < p' 2 . (7) 

Note that if the potential V contains an exponential ip dependence, as in hybrid inflation in supergravity 0], condition 
(0) can be violated and the instanton can have different from Eq. (^J) asymptotic behavior. In what follows we assume 
potentials which respect Eq. (Q). In the opposite case it is difficult to draw conclusions which are independent of a 
concrete form of a potential V(<p). 

The above discussion is true for any tpo for the monotonously increasing potentials. If the potential has a local 
minimum (false vacuum) dV/dip changes sign and under certain conditions [|l5| can stop the motion of tp. If tp' returns 
to zero precisely at cr/ one obtains the Coleman-De Luccia jl4| solution. In that sense the CD solution is a special 
(non-singular) case of the HT instanton. 

Some equations would be convenient to write in terms of a conformal time r, related to a by dr = da/b. Close to 
the origin a — e T (l + 0(e 2T )), and a — > corresponds to r — > — oo. 

Close to the singularity (which we assume at r = r/) da/dr = C(cr/ — cr) 1 / 3 and (r/ — r) = ^(o - / — cr) 2 / 3 which 
gives 



/ 3 2C 
<p = <Pf-J — 1n(Tf-T), b = C^j — {r f -t) . (8) 
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III. EUCLIDEAN QUADRATIC ACTION 



We will use the results of a well developed theory of gauge invariant cosmological perturbations to obtain the 
quadratic part of the Euclidean action. In |l3| ], in particular, cosmological perturbations for the Lorentzian version 
of the theory defined by the action Eq. (|l|) for closed universe have been considered. 

Since vector and tensor perturbations decouple in our model, we will be interested only in scalar perturbations, 
which are parameterized as follows 

ip -> <p(r]) + Sip (9) 
ds 2 = a(r/) 2 [-(l + 2A)drj 2 + 2B\ l dx l dr ] + ((1 - 2^)7^ - E^dx'dx 1 } , 

where 7^ is a metric on the unit 3— sphere and A, B, E, tp and Sip are small quantities. 

After the expansion of the perturbations by harmonics on the unit 3— sphere the Eq. (B18) from jl3| implies that 
the quadratic part of the action is represented as 

4 2) = £ {p2 -t KW I [(^) 2 + + ^(-u - p 2 + m<g\ *, , do) 

, Z J L UTj Z (f^ J 

where „ = d/dr/, K = +1 for a closed universe, p = 1,2, ... is an integer labeling different spherical harmonics and q 
is a gauge invariant quantity constructed from the perturbations^ 

q=—[iP-^(B-E >r ,)}. (11) 
K(p iT] a 

Now we analytically continue the Lorentzian quadratic action Eq. (|l0| ) to the Euclidean region taking rj = —it. We 

(2) (2) 

have S L = iS E where quadratic part of the Euclidean action is represented as a sum of following terms 
g(2) _ {p 2 - 4g)7r 2 [\ ( dq^ 2 ( K j 1 ^ 2 ^ 



S E 



+ M +M-V+P dr, (12) 

dr 2 (p 



where ' = d/dr. We see that quadratic action for the p = 2 mode vanishes. Another problem is that for the lowest p = 1 
(homogeneous) harmonic the overall sign of the quadratic action is "wrong" . We assume that analytic continuation 
q — > — if? is performed while integrating over this mode (compare fl6| , |ll1,|l2"f ) . 

The equation for the mode functions, which diagonalize the action (|12|), has form of the Schrodinger equation 

-^2l + U 1^ E 1^ ( 13 ) 

with a potential U depending on the background fields 

U^^ 2 + ^(-y+(p 2 -AK). (14) 
2 ip 

With the help of the equations of motion for the background quantities this potential can be written in the form 

b 2 2b 4 ,dV, 2 , 2 . 8bbdV l2 d 2 V 2k, 2 k 2 

17 = v + + (p - 4) - - 5 V - ^ bv - ^ ■ (15) 

One can check that close to the origin a — > where r — > —00 the first four terms contribute to the leading behavior 
and the potential tends to a positive constant 

U=p 2 + 0(c 2T ) . (16) 
Correspondingly the regular branch of the wave function dies out exponentially 



1 Since different harmonics decouple in what follows we shall omit the index p. 

2 We denote the scale factor a in the Euclidean region as b and assume K = +1 in what follows. 
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q = 7V e V^ r (l + 0(e 2T )) , (17) 

where N is a normalization constant. Note that E < p 2 corresponds to bound states and E > p 2 to scattering states. 

Towards the singular instanton end (as r — > rf) for sufficiently slowly growing scalar field potentials (condition (|t])) 
only the first and the last terms dominate in Eq. ( |l5| ) and give 

U= W ^+0((r f -rf). (18) 

If the potential V grows faster, U can change its asymptotic behavior and even sign, which might lead to catastrophic 
quantum mechanical behavior. 

The regular branch of the wave function in the potential Eq. ( |l8|) close to tj behaves as 

qnijt-Tf/*. (19) 

Note that in the derivation of the quadratic action Eq. ( |l2| ) we neglected the surface terms, which appear after 
integration by parts. The behavior Eq. ( |l7j ) and Eq. ( |l9| ) justify this procedure. 

In the intermediate region —oo < t < t/ the potential U can be positive as well as negative depending on the 
details of the scalar field potential V(ip), the initial value tpo and the harmonic number p. 

After rescaling q = f/cp and integrating by parts we obtain the action 



— 



which is obviously non-negative for p > 1. For p = 1 we have to investigate the problem numerically. 

IV. NUMERICAL RESULTS 

Let me describe our numerical results. The main idea is to study the zero energy wave function. Nodes of zero 
energy wave functions usually mean that there are bound states with negative energies. The number of nodes counts 
the number of bound states. 

For us it was convenient to solve the equations in ex. The a analogue of the Eq. ( |l3| ) has the form 

d 2 V d rrr7 E, 

^ = -6^ +(W -6^' (21) 

where 

Tir b 12 2 ,dV \ 2 p 2 -4 8b' dV d 2 V 2k k , 2 

w = A v + ^ ) + — w-^ y --^ ■ (22) 

One finds that the regular at a — > branch of the wave function behaves as 



q = N<7VP^(l + 0(a 2 )) , (23) 

where N is a normalization factor. In what follows we shall use k = 1 units. 

Starting at a = e with e = 10~ 3 4- 10 -4 we integrated simultaneously the background equations (||) with initial 
conditions (0) and the Schrodinger equation (^l|) with initial condition (|2^) for p = 1, E = 0, with different scalar field 
potentials and different initial values tpo. We investigated quadratic and quartic (monotonous) scalar field potentials 
V (up to parameter values V w Mpi) and a potential with the false vacuum. The shape of the potential U in the 
intermediate region depends on the initial conditions and details of the background scalar field potential V(tp). 

Typical results for the monotonous potentials are shown in Fig. |l| and Fig. for the case of a quartic scalar field 
potential V(<p) = ±A(1 + </). For some initial data the potential U is negative in the intermediate region, but in 



3 This fact was independently noticed in Q. 
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all cases the zero energy wave function has no nodes, i.e. there are no bound states with negative energies for these 
potentials. 

Next we investigated a potential 

2 

777 

V = —(p 2 (p-vf+Bp i ) (24) 

with m 2 = 2, v = 0.5, and B = 0.12 (compare ||). This potential has a local maximum at ipt op — 0.31250 and a 
local minimum (false vacuum) at ipf — 0.3571429. One finds the CD solution for tpo = p* = 0.1123575. It was shown 
in |Tl] ] that there are no negative modes about the CD solution. If one starts integration with ip < ip* one gets 
an "overshoot" singular solutions [HT instantons). We found that for these parameter values HT instantons have a 
negative mode as is shown in Fig. [| and Fig. ^[ For the ipo > <pf the potential is monotonous and no negative modes 
are found about corresponding HT instantons. For ip+ < ip Q < <p t one finds the singular solutions or Bousso-Lindc 
double- bubble instantons ||. 



V. CONCLUDING REMARKS 



To conclude, we demonstrated analytically that the quadratic action for the HT instanton is non-negative for the 
harmonics with p > 1. For the lowest p = 1 (homogeneous) harmonic we do not have a general proof. We investigated 
corresponding Schrodinger equation numerically for different scalar field potentials and different scalar field initial 
values and in some cases found a bound state. We argue that there is a negative mode about the HT instanton for 
certain potentials and certain initial conditions. 

In H it was suggested that the HT solution be considered as a constrained instanton. It is an open question how 
one can incorporate this idea into the stability consideration. 

As a starting point of our analysis the expression for the quadratic action has been taken from the work [ p"3"[ , 
where it was obtained via Hamiltonian reduction. It is interesting to note that using the Lagrangian formalism one 
obtains only a sign indefinite quadratic action [jl7|,[llj , whereas the Hamiltonian reduction allows one [fy]jl3) to derive 
a quadratic action which has negative overall sign only for the lowest (homogeneous) harmonic. This problem does 
not arises for the K = universe. 

Let me mention that the quadratic action in the form Eq. (|l^) cannot be straightforwardly used to study the 
double- bubble instantons of Bousso and Linde [^j, since these solutions have points ip = in the intermediate region. 

These questions need further investigation. 
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FIG. 1. 95(f) and b(a) for the HT instanton for the potential V = |A(1 + ip 4 ) with A = 0.01 and ipo = 0.1. 




FIG. 2. The potential U and the zero energy wave function q (for p — 1 harmonic) corresponding to the instanton shown 
in Fig. ^. The wave function has no nodes. 
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FIG. 3. <fi(a) and b(a) for the HT instanton for the potential V = ^-(^ 2 (v - vf + Bip 4 ) with m 2 = 2,u = 0.5, B = 0.12, 
and ipo = 0.11 < 93*. 
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FIG. 4. The potential U and the zero energy wave function q (for p = 1 harmonic) corresponding to the HT instanton 
shown in Fig. ta. The wave function has a node. 
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